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The general extreme limit of the double-Reissner-Nordstrom solution is worked out in explicit 
analytical form involving prolate spheroidal coordinates. We name it the combined Majumdar- 
Papapetrou-Bonnor field to underline the fact that it contains as particular cases the two-body 
specialization of the well-known Majumdar-Papapetrou solution and Bonnor's three-parameter 
electrostatic field. To the latter we give a precise physical interpretation as describing a pair of 
non-rotating extremal black holes with unequal masses and unequal opposite charges kept apart by 
a strut, the absolute values of charges exceeding the respective (positive) values of masses. 
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I. INTRODUCTION 

The Majumdar-Papapetrou (MP) [1, 2] electrostatic solution of the Einstein-Maxwell equations can describe in 
the axisymmetric case a system of aligned non-rotating charged masses which are in neutral equilibrium due to 
balance of the gravitational and electric forces. Although it was first thought that this solution "is not of great 
physical interest" [3] , later on the well-known paper by Hartle and Hawking [4] , in which the charged sources of the 

■ MP spacetime located on the symmetry axis were shown to be the extreme Reissner-Nordstrom (RN) black holes 
' [5, 6], quickly converted the MP solution into an important ingredient of the modern black-hole physics, with some 

distinctive thermodynamical properties [7]. Since the MP black-hole spacetime formally belongs to the Weyl family 
of electrostatic solutions known from 1917 [8], it is quite surprising that it was only discovered three decades after 

, Weyl's work, besides, first in the framework of non-symmetric, conformastatic spacetimes [9]. This could probably 
^ ' be explained by the fact that already in the two-body case the performing of the extreme limit in the Weyl solution, 

| needed for obtaining the corresponding MP metric, is technically so tangled that it has never been attempted in the 

■ literature. In this respect it would certainly be of interest to examine the limit of extreme black holes in the five- 
parameter double-Reissner-Nordstrom (DRN) solution [10, 11] which describes the non-linear superposition of two 
arbitrary RN fields and hence contains the two-body electrostatic Weyl spacetime as a particular case. By considering 
such a limit one could achieve the following two goals: (i) to show how the MP solution for two extreme RN black 
holes naturally arises from the general DRN spacetime, and (ii) to identify and analyze the remaining extreme cases. 

In the present paper, following the main ideas of Rcf. [12] dealing with a pair of stationary electrovac extreme 
sources, we will work out the general extreme limit of the DRN solution by making use of the canonical set of 
parameters {a n ,/3i} in terms of which earlier in the paper [10] were written both the general metric for N aligned 
RN sources and its N = 2 specialization representing two black holes. By the term 'extreme', with regard to a 
non-isolated RN black hole, we shall mean, as usual, the marginal black-hole state previous to the formation of a 
naked singularity. Curiously, it turns out that the physical representation of the DRN solution obtained in [11] by 
introducing the Komar masses and charges [13] as arbitrary parameters, is inappropriate for achieving such a goal 
since it produces hardly manageable expressions while performing the limit o\ — > 0, 02 — > 0. However, formulae 
relating various mathematically equivalent parameter sets that were derived in the paper [11] will be used by us for 
clarifying the physical meaning of the resulting extreme solutions. 

Our paper is organized as follows. In Sect. II we first provide the reader with a necessary information about the 
canonical parametrization of the DRN solution and then elaborate the limit of two extreme black holes, writing it 
down in terms of two determinants of fifth order. In Sect. Ill we consider expansions of these determinants, arriving as 
a result at the combined Majumdar-Papapetrou-Bonnor spacetime which contains as special cases the MP two-body 
solution [1, 2] and Bonnor's three-parameter electrostatic solution [14]. Here we also show that the Bonnor solution 
is just the electrostatic analog of the well-known Kerr-NUT spacetime [15], thus obtainable from the latter via a 
complex continuation of the parameters [16], and we interpret it, in the line of papers [17, 18], as describing two 
extreme black holes with unequal masses and unequal charges (the absolute values of charges being greater than the 
respective positive values of masses) which are separated by a strut. Discussion and conclusions are left for Sect. IV. 
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II. THE EXTREME DRN SOLUTION IN TERMS OF DETERMINANTS 



We recall that the general DRN solution arises from the axis data (short for values of the Ernst potentials [19] on 
the upper part of the symmetry axis) of the form [10] 
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where the parameters e;, 0i and I = 1,2, can take on arbitrary real values or occur in complex conjugate pairs 
ei = e2, 0x = 02, fi = ]i (all of them simultaneously) due to the reality of the potentials £ and <3> in the electrostatic 
case. Note that these parameters define only five physical real quantities because the axis expressions of the Ernst 
potentials (1) include an arbitrary shift along the symmetry z-axis which can always be fixed at some real value. The 
total mass M to t and total charge Qtot of the system are given by the formulas 



M tot = -(ei + e 2 )/2, Q tot = fx + f 2 . 



(2) 



The application of Sibgatullin's integral method [20, 21] to the axis data (1) yields the following expressions for £ 
and $ in the whole space [10]: 



E+/E-, 
1 

±1 



$ = F/E. 



E± 



1 

n 



l 

»*4 



±1 






a x - 0x 
ri 

ai - 2 
hijai) 
ax - fix 
h 2 (ax) 



r n = 



ax - f>2 
Vp 2 + {z- a n ) 2 , 



a 4 - fix 

hx{dA) 
ai - fix 

ai - fj 2 
hi(a n ) -- 



F 






ei + 2f{f(a n ) 





■ ■ ■ f{ati) 


n 


ri 


ax - fix 


a A - fix 




ri 


ax - fj 2 


cti - fJ 2 


hx(ax) 


hx(a 4 ) 


ax - fix 


ai- fix 


h 2 {ax) 


h 2 (a 4 ) 


ax - fj 2 


OLi- fj 2 


f{0tn) 


2 , 
\p Jl 



0i 



(3) 



where instead of the initial set of parameters {e/, /?;,//}, a mathematically equivalent parameter set {a„,/3;} is used, 
the new constants a ni n — 1,2,3,4, taking arbitrary real values or occurring in complex conjugate pairs. The form 
of the 'old' parameters ei and /; in terms of a n and fii is the following [10]: 
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The constants a n define the location of sources on the symmetry axis. When all a n are real, these define two black 
holes; the complex conjugate pairs of as determine hyperextreme constituents. Therefore, it is clear that for passing 
from the case of two subextreme black holes to the case of extreme black-hole constituents we should work with the 
real- valued a n , assigning them first of all some order; the order adopted here is 



ax < a 2 < a% < ai. 



(5) 



Then the parts p — 0, a 2 < z < ax and p — 0, ai < z < a 3 of the symmetry axis will define the horizons of the upper 
and lower black holes, respectively (see Fig. la). The limiting procedure consists in tending a 2 to ax, and a 4 to a% 
in formulas (3) and (4). The limit a 2 = ot\, a 4 = a% in (4) does not exhibit any difficulty; however, it is necessary 
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to be careful while evaluating /i and / 2 in order not to omit any sign that these quantities might take. Then, in the 
extreme limit we readily get from (4) 
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the corresponding expressions for e\ and e 2 being 
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As will be seen later, the concrete choice of ei and e 2 in (6)-(7) determines whether the charges of black holes have 
equal or opposite signs. 

The limiting procedure in the expressions defining Ernst potentials (3) requires the application of l'Hopital's rule 
because, with a 2 = a\ and 0:4 = a 3 , the determinants E± and F vanish. However, we can write down immediately the 
final result for £ and $ if we observe that using of l'Hopital's rule in (3) is equivalent to the application of Sibgatullin's 
method to the axis data (1) under the supposition that the resulting solution will involve a pair of as of multiplicity 
two, something that was already done in the paper [12] in a more general context of stationary electrovac spacetimes. 
Therefore, using the results of [12] and taking into account the reality of the Ernst potentials in the electrostatic case, 
below we can write down the final form of £ and $ in the required extremal limit: 
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where, according to the conventions of paper [12], partial differentiation in the determinant F and in the matrix A is 
applied only to the quantities f(a n ), (a n — /3;) and functions r„, the constants e/, Pi, /; being assumed independent 
of ai and a 3 ; thus, for instance, 
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The form of the corresponding metric coefficients / and 7 entering Weyl's line element 

As 2 - rV 7 (V + dz 2 ) + pW] - fdt 2 , 
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is obtainable from formula (14) of [12], yielding 
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Therefore, we have found a representation of the extreme limit of the DRN solution in terms of 5 x 5 determinants. 
It should be pointed out that once the general expressions for the Ernst potentials and for the whole metric are 
obtained, these can be further reparametrizcd by employing any of the equivalent parameter sets that one might have 
at his disposal. At present moment we go on using the set {a„,ft} which, after the limiting procedure has been 
carried out, reduces to only four arbitrary constants, ai, a 3 , ft, ft, representing three physical quantities (recall 
that the sum of ai and a 3 can be assigned any real constant value). The form of fi and ei which enter the constant 
objects hi(a n ) is defined in terms of a n and ft by formulas (6) and (7). Our next major objective will be rewriting 
the obtained general expressions for the potentials £, $ and metric functions /, 7 in a simpler form by substituting 
formulas (6), (7) into (8), (11) and expanding the determinants. 



III. THE GENERAL EXTREME SOLUTION IN SPHEROIDAL COORDINATES 



Although the procedure of the expansion of determinants in formulas (8) and (11) involving the general values of fi 
and ei might look complicated at first glance, this is not really the case because the corresponding analytical computer 
processing does not exhibit much difficulty and leads straightforwardly to the desired result. We only mention that it 
is advantageous to choose the origin of coordinates in such a way that ai + a 3 = (see Fig. lb), and also introduce 
prolate spheroidal coordinates (x, y) via the formulas 



ai=a, a 3 = -a, x = ^-{r 3 + n), y = ^-(r 3 -ri), 

la la 

where a is a positive real constant. The form of n, r 3 and z in terms of x and y is 

n=a(x-y), r 3 =a(x + y), z = axy. 
In the new coordinates the line element (10) rewrites as 
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Below we give the final expressions for £, $, /, 7 at which we have eventually arrived after the substitution of (6), 
(7) into (8), (11) and simplifications: 
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Formulas (15) are fully equivalent to formulas (6), (7), (8), (11) of the previous section, but due to their concise 
form they are of course by far more advantageous as a representation of the general extreme DRN solution than the 
determinantal expressions. In our opinion, the most proper name for the field defined by (15) would be the combined 
Majumdar-Papapetrou-Bonnor solution, and in what follows we will justify this name by considering two particular 
specializations of the solution (15) arising from two different relations between the quantities e\ and e 2 . 



A. The Majumdar— Papapetrou two— body solution 

When ei = e 2 = e, e = ±1, we obtain from (15) 
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and these Ernst potentials define the Majumdar-Papapetrou solution for two extreme Reissner-Nordstrdm black 
holes. To see this better, let us also write out the corresponding metric functions / and 7: 

f= a^l t e2 , = 1> (17) 

whence it follows that the solution (16) is conformastatic [9]. Moreover, in the particular case under consideration 
one finds from (2), (6) and (7) that 

ei + e 2 = -2e(/i + / 2 ) <f=> Mi + M 2 = e(Q 1 + Q 2 ), (18) 

Mi being the Komar masses and Qi the Komar charges of black holes. To show that this implies the equalities 
Qi = eMi characterizing the Majumdar-Papapetrou solution, one has to make use of the quantities o\ and o~ 2 of 
Rcf. [11]. Then, expressing the separation R from the equation o\ = and substituting the result into the equation 
a 2 = 0, one can easily solve the latter for M 2 ; this M 2 being substituted into (18) finally gives Qi = eM\. In 
analogy one obtains Q 2 = eM 2 . Therefore, the equivalence of the solution (16) to the N — 2 specialization of the 
Majumdar-Papapetrou spacetime is evident. 

The results of the paper [11] also permit one to relate the parameters /3i and f3 2 entering (16) to the Komar masses 
and to the distance a = R/2 of each extreme source from the origin of coordinates: 

A + P2 = -Mi - M 2 , A/3 2 = -a(a + M x - M 2 ) (19) 

(the subscripts 1 and 2 label, respectively, the lower and upper constituents), so that the polynomials E± and F 
defining a representation of the two-body Majumdar-Papapetrou solution in prolate spheroidal coordinates can be 
rewritten as 

E± - a(x 2 - y 2 ) T (Mi + M 2 )x ± (Mi - M 2 )y, 
F = e[(M 1 +M 2 )x-(M 1 -M 2 )y], (20) 

thus allowing for a slight simplification of formulae (16) and (17). 

It is worth noting that, as it follows from a comprehensive analysis of MP spacetimes carried out by Hartle and 
Hawking [4], Eqs. (16), (17) describe a regular geometry of two charged extreme black holes only in the case Mi > 0. 
If one of the masses is a negative quantity, the solution develops naked singularities off the symmetry axis, even when 
the total mass of the system is positive. The above said is well illustrated by Fig. 2 in which we have plotted singular 
surfaces for the following two particular parameter choices: (i) Mi = 2, M 2 = —1, a = 3, and (ii) Mi = —2, M 2 = —1, 
a = 3. For that reason, throughout the present paper only the constituents with positive Komar masses are referred 
to as black holes. 

B. Bonnor's solution: the electrostatic analog of Kerr— NUT spacetime 

The other relation between ei and e 2 is ei = — e 2 = e, e = ±1, and in this case we get from (15) 

E+ &2l _ N± 

° 77i ' XT, ' J 7712 ' D _ 



' J EV " a 8 (/3i -p 2 f{x 2 -y 2 ) 4 ' 



6 



E± = KA-ftjx-^-A^i^-a 2 )] 

x [a{fa - fa)x + (a 2 ~ fafa)y ± (a 2 - /3 2 2 )], 

F = a(2a 2 -fa 2 - /3 2 )x - {fa + fa)(a 2 - fafa)y, 

N = a 2 (fa - fafx 2 - (a 2 - fafafy 2 + (a 2 - (3 2 )(a 2 - fa 2 ). (21) 

This is the three-parameter electrostatic solution first obtained by Bonnor [14], though in a different form and not 
revealing that it was the electrostatic analog of the well-known Kerr-NUT spacetime [15], i.e., obtainable from the 
latter via Bonnor's procedure of a complex continuation of the parameters [16]. In the particular case 2a 2 = fa 2 + fa 
it further reduces to the solution for a massive electric/magnetic dipole [22]. The relation between our parameters a, 
fa, fa and the constants a, A, B, C of the original Bonnor's paper [14] is given by the formulae 

a = a, fa=a(A + B), fa = -a(A + C), BC = A 2 - 1. (22) 

The simplest way to show how (21) can be constructed from the Kerr-NUT solution is to use the expression of the 
corresponding Ernst potential from Ref. [23] (formulae (3) of [23] with q = b = 0) and the Bonnor representation of 
the axisymmetric electrostatic problem (equations (2.9) and (2.10) of [14]). Then, after changing a — s> ia, v — > iv in 
the Ernst potential £ of the Kerr-NUT solution, and also in its complex conjugate expression £, one will arrive at 
two real potentials X and Y whose product will give precisely the electrostatic Ernst potential £ from (21), while 
the difference of X and Y will give the doubled value of the potential $ from (21). The relation of the parameters 
a, fa, fa entering (21) to the Kerr-NUT parameters m, a, v (the latter two after their complex continuation) is the 
following: 

a = \/m 2 + a 2 — v 2 , fa = —m + a + v, fa = —m — a — v. (23) 

It should be pointed out that solution (21) has important aspects of its physical interpretation that need to be 
clarified. In the original paper [14] it was only referred to as describing a mass endowed with both electric charge 
and dipole moment. At the same time, during the last decade there appeared some new interesting results shedding 
additional light on the interpretation of Bonnor's solutions. Thus, for instance, in Ref. [17] the electric/magnetic 
dipole spacetime [22] was shown by Emparan to arise from two charged extreme black holes with equal masses and 
equal but opposite charges, and later on this result was generalized by Liang and Teo [18] to the case of the non-zero 
net charge of black-hole constituents. However, whereas in the paper [17] a correct relation between the charges and 
masses of extreme black holes was established, the authors of [18] presented erroneous relations between the masses 
and charges of the constituents 1 and, besides, they did not identify the electrostatic analog of the Kerr-NUT metric as 
the Bonnor three-parameter solution [14] (but, importantly, they correctly indicated, on the one hand, the inequality 
of masses, and, on the other hand, the inequality of charges in their two-body configurations). A crucial mistake 
made by Liang and Teo was to assume that the extreme black holes in their Einstein-Maxwell dihole solution were 
of the MP type, i.e., verifying Q 2 — M 2 , with the only distinction that the charges had to have opposite signs. 

The correct formulae relating the charges and masses of the extreme black-hole constituents of the solution (21) 
are readily obtainable by equating the quantities o\ and a 2 of the paper [11] to zero and excluding the MP case; we 
thus get 
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(R + M 2 ) 2 -M 2 {R + Mi) 2 - M 2 

Ql - eMl \l (R-M 2 ) 2 -Mr ^--^(R-Mtf-M*' (24) 

whence the Liang-Teo conjectured relations follow only in the limit R — >• oo of infinite separation. One can also see 
that Emparan's result for the Bonnor dipole solution is recovered from (24) when Mi = M 2 . 

Formulae (24) clearly demonstrate that the charges of extreme constituents in the solution (21) always have opposite 
signs, provided that Mj > 0. Another important conclusion following from (24) is that the absolute values of these 
charges are greater than the respective values of the masses, \Qi\ > Mj. This curious property of the non-isolated 
extreme charged black holes was first established by Emparan in the case of identical (up to the sign of the charge) 
black-hole constituents [17]. The present paper extends Emparan's discovery to the case of extreme constituents of 
the non-MP type with unequal positive masses. Mention here that the negative masses of the constituents are likely 
to be excluded on the same grounds as in the MP case - they are accompanied by naked singularities, and in Fig. 3 



1 For example, a correct expression for the masses of extreme black holes in the absence of dilatonic field is m ± {l/a)\/m 2 + a 2 — I 2 (in 
notations of [18]), and not mil as was given by Liang and Teo. 
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we have plotted singular surfaces for the same choices of the parameters Mi, a which were earlier used for plotting 
Fig. 2. 

We would like to conclude this section by noting that the physical quantities are related to the parameters a, fix, 
02 of Bonnor's solution by the formulae 

(a + /3i)(a + /3 2 )(/3i+/3 2 -2a) 
IVIi = 
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while the solution (21) can be rewritten in terms of the physical parameters R, Mi, M 2 by means of the substitutions 

a = ±R, pi = -^{D + Mi + M 2 ), /3 2 = ^(D-M 1 -M 2 ), 



= (R-M 1 -M 2 )[(R + M 1 )^-Mj] 

V R-Mi + M 2 ' { ' 

The above relations between the parameter sets are a direct consequence of the results of Ref. [11]. 



IV. DISCUSSION AND CONCLUSIONS 



In the present paper we have elaborated the general extreme limit of two non-rotating charged black holes starting 
from the five-parameter DRN solution in the canonical representation. This limit is represented by the combined 
Majumdar-Papapetrou-Bonnor field whose particular specializations are the two-body MP spacetime written in 
prolate spheroidal coordinates (when the charges of black holes have equal signs) and Bonnor's electrostatic solution 
(the case of oppositely charged masses). The extreme limit involves in general three arbitrary real parameters which 
can be associated with the Komar masses of the black holes and the separation distance. At the same time, whereas 
the extreme constituents of the MP solution are in equilibrium which is distance-independent, the equilibrium states 
within the Bonnor solution do not exist, as was shown by Liang and Teo [18], so the constituents in this particular 
solution are kept apart from each other by a supporting strut [24] . An important characteristic feature of the extreme 
black holes in the Bonnor solution demonstrated by us is the inequality \Q%\ > Mi between the individual charges and 
masses, and it generalizes Emparan's result previously obtained for a particular case of identical constituents [17]. 
This means, taking into account the equality Q 2 = M 2 valid for extreme constituents of the MP solution, that in the 
two-body systems of charged extreme black holes the absolute value of the individual Komar charge can never be less 
than the respective individual Komar mass. 

It is worth noting that formulas (16)-(20) describing the two-black-hole MP spacetime in prolate spheroidal co- 
ordinates are well suited for the illustration that the extreme black-hole constituents are not the only axisymmetric 
objects in equilibrium covered by the MP class of solutions. Indeed, while working with the x and y coordinates, it 
is quite customary to explore both the prolate and oblate variants of the coordinate system, so in (16)-(20) we can 
pass to the latter (oblate) variant by the transformation 

x — > ix, a — > — ia, y — > y, (27) 

which, however, should not affect the reality of the Ernst potentials £ and $. From (20) it is easy to see that the 
condition for £ and $ to continue to stay real under (27) is the equality of the masses Mi and M 2 . Therefore, 
after setting Mi = M 2 — \M and passing in the MP solution for two extreme black holes to the oblate spheroidal 
coordinates by means of (27) we arrive at the following transformed metric and electric $ potential: 



d S 2 = a 2 /- 1 



(■'■ 2 + /r)(^ I + I ^]+^ + i)(i-. v W 



-fdt 2 



= a 2 (x 2 +y 2 ) 2 = eMx 

3 ' [a{x 2 +y 2 ) + Mx] 21 a(x 2 + y 2 ) + Mx K ' 

The above formulas represent what we would call the simplest MP ring, whose mass and charge are both equal to M 
in absolute value, and in the paper [25] this solution arose as a particular static limit of a more general stationary 
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electrovac spacetime. For us, the most interesting point about the solution (28) is that although its mass and charge 
satisfy, like in the case of an isolated extreme RN black hole, the relation Q 2 = M 2 , it is not already spherically 
symmetric because it possesses higher multipole moments [25]. 

If we now return to Weyl's cylindrical coordinates and superpose the ring solution with the extreme RN field, we 
shall arrive at a particular N = 3 specialization of the MP spacetime of the form 

ds 2 = (l + ^+^) 2 (dp 2 +dz 2 +pW)-(l + ^ + ^)- 2 dt 2 , <&= £{qpr + qph) 



1+lpr+lph' 

M ( 1 1 \ TO 

A ~ ~2 \^p^ + {z + iaf + VP 2 + (* - iaf J ' V ' P 1 + (z - b) 2 ' 

where m is the mass of the black hole and b its position on the symmetry axis (the subindices r and h refer to the 
ring and the black hole, respectively). It is remarkable that the ring and and the extreme black hole are forming 
a distance-independent equilibrium configuration which is similar to the one of the two-black-hole MP solution. 
Moreover, a trivial modification of the potential tjjh in (29) leads to the equilibrium states between two MP rings 
that are also distance-independent! Therefore, actually we have at hand three different geometries representing the 
two-body equilibrium configurations independent of separation and involving the equality Qi — eMj. In our opinion, 
this might be regarded as a possible indication that the classical equilibrium condition MiM 2 = Q1Q2 for charged 
point particles must also hold in general relativity if the effects of the induced electric field [26, 27] are not taken into 
account. 
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FIG. 1: Location of the black-hole constituents on the symmetry axis: (a) the subextreme case; (b) the case of two extreme 
constituents. 



10 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 

P P 

FIG. 3: Formation of the singular spheroidal surfaces in the case of the Bonnor solution when one or two constituents have 
negative masses. 



